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Abstract. — The temperature and magnetic field dependence of the specific heat Cp(T,H) in the 
superconducting mixed state as well as the upper critical field Hc2{T) have been measured for poly- 
crystalline Y2;Lui_a;Ni2B2C and Y(Nii_j;Pty)2B2C samples. The linear-in-T electronic specific heat 
contribution 'y{H) ■ T exhibits significant deviations from the usual linear-in-_ff law for all x and y the 
transition metal site (T) resulting in a disorder dependent negative curvature of 7{-ff). The deviations 
from that linear behaviour of our unsubstituted samples are the largest reported so far for any super- 
conductor. The Hc2{T) data point to the quasi-clean limit for (Y,Lu)-substitutions and to a transition 
to the quasi-dirty limit for (Ni,Pt)-substitutions. The -f{H) dependence is discussed in the unitary 
d-wave as well as in the quasi-clean s-wave limits. From a consideration of j{H) data only, d-wave 
pairing cannot be ruled out. 



INTRODUCTION. - The rare earth (R) transition metal (T) borocarbide family 
(i?C)„T2B2 (i? = Y, Lu, So, Th, La; T = Ni, Pd, Pt; n = 1 or 2) contains superconductors 
with relatively high transition temperatures Tc up to 23 K The coexistence of super- 

conductivity and magnetism for members of this family where R are magnetic rare earth ions 
such as Dy, Ho, Er, ... , has stimulated numerous studies of their thermodynamic and trans- 
port properties in the superconducting as well as in the normal state. At first glance, most of 
those results support a classification of these materials as intermetallic phonon mediated su- 
perconductors with a moderately strong coupling strength. However, clean i?Ni2B2C samples 
exhibit also some features unexpected for ordinary s-wave superconductors. We emphasize 
the unusual shape and the strong disorder dependence of the upper critical field Hc2 (T) and a 
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nearly scaling of the electronic specific heat Ces (T) in the superconducting state compared 
with exponential behaviour for ordinary s-wave superconductors . 

According to Nohara et al. Q the isoelectronic T-substitution does affect strongly the field 
dependence of the linear-in-T electronic specific heat contribution ^{H) ■ T in the mixed state. 
Thus, for an Y(Nio.8Pto. 2)2620 single crystal 7(i?) oc H has been found, while a square-root 
law was observed for a pure YNi2B2C single crystal and for polycrystalline LuNi2B2C [|j 

7(7?) /7iv (X V^/^c2(0), (1) 



where is the Sommerfeld constant in the normal state. Although the observed 7(i?) c>c v H- 
law for YNi2B2C and LuNi2B2C was regarded initially as evidence for d-wave pairing |^,^, 
the disorder related transition from a ^/H to a lincar-in-_ff dependence was subsequently used 
to rule out d-wave superconductivity in non-magnetic borocarbides However, to the best 
of our knowledge, systematic investigations of this problem in a broader concentration range 
for Y(Nii_j,Pty)2B2C are lacking. Since isoelectronic substitutions in the RC charge reservoir 
are expected to produce much weaker disorder than those in the TB network we studied also 
the closely related Ya;Lui_a;Ni2B2C system for the sake of comparison Q]. By changing both 
compositions, x and y, deeper insight should be gained on how does the disorder affect the 
field dependence of the specific heat Cp(r, iJ), the shape and the magnitude of Hc2{T), as well 
as the nature of the pairing state. 

EXPERIMENTAL DETAILS. - Polycrystalline Y:rLui_^Ni2B2C with x = 0, 0.25, 
0.5, 0.75, 1 and Y(Nii_j,Ptj,)2B2C samples with y = 0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.5 and 0.75 
were prepared by a standard arc melting technique. Powders of the elements were weighted 
in the stoichiometric compositions with a surplus of 10 wt.% boron to compensate losses of 
boron during arc melting. The powder was pressed to pellets which were melted in argon 
gas on a water-cooled copper plate in an arc furnace. To get homogeneous samples, they 
were turned over and melted again four times. After the melting procedure the solidified 
samples were homogenised at 1100 °C for ten days. The specific heat was measured in the 
range 4.2 K < T < 20 K and for magnetic fields ^qH < 8 T using a quasi-adiabatic step 
heating technique. The heating pulses were generated by a strain gauge heater and T was 
measured with a Au-Ge thin film resistor. The temperature dependence of the upper critical 
field Hc2{T) was determined by taking Tc from the jump of Cp in the particular field. 

RESULTS AND DISCUSSION. - To illustrate typical specific heat behaviour, the 
Cp/T vs data at H = oi the Yj;Lui_a;Ni2B2C series and the corresponding curves for 
HqH < 8 T of the pure Y sample {x = 1) are shown in fig. 0. Measurements at 8 T were used to 
analyse the normal state specific heat Cp = jnT + PdT^, where 7jv is the Sommerfeld constant 
and PdT^ is the Debye contribution. The 7jv was determined by extrapolating the Cp/T vs 

curves of the high field data in the normal state to T — > 0. In this way we obtained 7Ar = 
20.4 {x = 0), 19.0 (x = 0.25), 18.3 (x = 0.5), 18.0 {x = 0.75) and 20.2 mJ/molK^ {x = 1) for 
our Y2,Lui_a;Ni2B2C series in good agreement with the data reported by several groups |p|-jl^ 
and 7Ar = 20.2 (y = 0), 20.2 {y = 0.05), 18.4 (y = 0.1), 16.4 {y = 0.15), 16.2 (y = 0.2), 16.9 
(y = 0.25), 15.3 (y = 0.5) and 15.0 mJ/molK^ (y = 0.75) for the Y(Nii_2,Pty)2B2C series. To 
determine j{H), the Cp/T vs curves in the particular field have been extrapolated in the 
same way to T — > from the data in the range 4.2 K < T < 7 K or up to the onset of the 
transition to the normal state (see fig. 

For all samples j{H) is a sublinear function of H. At first generalizing eq. (|^), the data 
were analysed by the expression 

7(iJ)/7^ = [H/H^^m'-^, (2) 
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where the fitting parameter (3 measures the subhnearity {i.e. a negative curvature) of 7(iJ) 
and Hc2{0) is the field where -/{H) reaches 7jv. We obtained /? = 0.66, 0.42, 0.42, 0.41 and 
0.56 ongoing from x = to x = 1 for Y^Lui_::,Ni2B2C and /3 = 0.56, 0.36, 0.35, 0.26, 0.18, 
0.20, 0.49 and 0.60 ongoing from y — Q toy = 0.75 for Y(Nii_^Pt„)2B2C, with uncertainties of 
^Pl P < 10% due to the smaU measured values 70 = ^{H = 0) ||l^ and due to the procedure 
used to determine jiH), as mentioned above. The dependence of P{x) is shown in the inset 
of fig. H and in fig. |j (left panel). /3 reaches the largest values for the bordering cases x — 
and 1 and becomes markedly smaller in between. We note that our curvatures for LuNi2B2C 
and YNi2B2C exceed slightly the value of /3 = 0.5 suggested by eq. (|l|) and that reported in 
Refs. mis). To the best of our knowledge the strong sublinearities for j{H), measured by the 
exponent /3, of the borocarbides under consideration are the largest reported so far for any 
superconductor. The P{y) behaviour is depicted in the inset of fig. ||and in fig. ^ (right panel). 
Large exponents P « 0.5 are observed for y = and y — 0.5, whereas /3 is significantly reduced 
for Pt-concentrations y inside this range and has a finite minimum value /3 = 0.18 at y = 0.2. 
This minimum value is at variance with the linear law for an Y(Nio.8Pto.2)2B2C single crystal 
reported in ref. [Q. We attribute that observation to a stronger disorder compared with our 
samples. Larger Pt-concentrations reveal even more pronounced (3 values, e.g. y — 0.75; 
P = 0.6. 

The observed j{H) oc H^^^-law with /3 « 0.5 raises the question whether an uncon- 
ventional pairing mechanism is responsible for this peculiarity since according to Ref. [ p6[ 
7(i?) (X VH is a signature for a nodal order parameter with d-wave symmetry (a somewhat 
larger value /3 = 0.59 has been found in Ref. |17[ |) while jiH) cx H describes an isotropic 
s-wave order parameter. According to Refs. m8, 19| Volovik's clean limit d-wave approach 
can be generalized to describe also strong impurity scattering. Then at low magnetic fields 
H <^ Hc2{0) the specific heat coefficient 'j{H) follows a. H\nH dependence: 



7(i?) = 70 -f JnD 



H 



In 



2^ 



gc2(0) 

H 



(3) 



where a and D are constants. Such a behaviour was observed for various disordered high- Tc 
cuprates and considered as evidence for d-wave superconductivity in the unitary scattering 
limit ||2^,^. At the same time its applicability to non-magnetic borocarbides under consider- 
ation was disclaimed (2^]. However, some of our data can be described equally well by eq. (||) 
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Fig. 2 — Magnetic field dependence of the specific heat contribution ^{H) of the vortex core electrons in the 
mixed state (H < Hc2) normalised by 7jv with //c2{0) as the upper critical field (see fig. W) for Ya;Lui_j;Ni2B2C 
(left panel) and Y(Nii_j,Pty)2B2C (right panel). The lines are fits according to eq. and the straight 
reference line corresponds to the usual linear-in-i? s-wave dirty limit behaviour. The insets show the curvature 
parameters fi{x) and /3(y) as defined in eq. (H). 



for H/Hc2{0) < 0.3 as well as by eq. (||) using intermediate values for /3 (0.2 to 0.35). This is 
shown in fig. ^ obviously, the H\nH behaviour is not very distinct from the power law at low 
fields fioH < 1.5 T. At higher fields the H\nH dependence may deviate since it was derived 
for low fields only For y = 0.25 and 0.5 we observed residual values 70 = 3.4 and 3.3 

mJ/molK^, respectively, which are substracted in fig. ^ The existence of a non-negligible 
7o is a feature predicted for a d-wave order parameter in the unitary limit p^ . Hence, d-wave 
pairing cannot be ruled out in non-magnetic borocarbides by considering ^{H) data only. 
While the deviation from the linearity of ^{H) is frequently ascribed to a shrinking of the 




Fig. 3 — Magnetic field dependence of 7(H) for Y(Nio.75Pto.25)2B2C. The solid line is a fit according to eq. 
(tSl). The dashed line is a fit according to eq. (H) with (3 = 0.20. 



vortex cores with magnetic field and to vortex core interactions j4[^,|2j], recent investiga- 
tions support the assumption of delocalized quasiparticle states around the vortex core to be 
responsible for this feature, similar as in d-wave superconductors ||2^]. However, there are sev- 
eral conventional, but anisotropic s-wave, superconductors which also exhibit deviations from 
the 7(if) (X H-lem in the clean limit, e.g. VgSi M, NbSes [| (/3 = 0.33), and CeRuj 0,1^. 
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Remarkably, a sublinear j{H) behaviour has been reported also for the novel "medium- Tc" 
superconductor MgB2 . In this general context recent ultrahigh-resolution photoemission 
spectroscopy measurements suggest that a highly anisotropic gap might be responsible for the 
mentioned above peculiarities in clean Ni borocarbides |29[ . By introducing disorder (y = 0.2) 
a complete isotropization of the gap was observed. Calculations of the density of states (DOS) 
at the Fermi level, N{0), in the mixed state with interacting vortices by Ichioka et al. JTzt 
revealed a — 0.33) dependence of 7(i?) for anisotropic s-wave superconductors. 

In this context it is noteworthy that our low-temperature Cp data on YNi2B2C (x = 1) 
point to a possible exponential behaviour of Ces{T) below 1/6 with A{0) /ksTc « 0.5. 
This supports an effective multiband s-wave picture like that proposed by Shulga et al. pO| . 
Noteworthy, a similar value for A(0) was observed in microwave measurements pl| . 

Like jiH), the upper critical field Hc2{T) can be described also by a simple scaling law 



Hc2{T) = H*2{'^-TlTc)^+", valid for 0.3 < T/T^. (4) 

Our values of the upper critical field Hc2{0) ~ 0.9i/*2 are reduced due to R-site substitution. 
A similar behaviour was found for the pronounced positive curvature of Hc2{T) near 
which is measured by the exponent a in eq. (Q), in contrast to the opposite statement of a 
nearly constant curvature The unusual positive curvature of Hc2{T) near observed 
here can be explained for superconductors in the clean limit by a significant dispersion of 
the Fermi velocities using e.g. an effective two-band model pO| ]. Tc and ^jy are reduced to a 
smaller extent, which has been ascribed to a slight reduction of the electron-phonon coupling 
constant A at intermediate x For Tc a dip near x = 0.7 is observed, in accordance 

with refs. (Tc ~ 14.6 K at a; = 0.75). The dirty limit region is not reached (which would 

be represented by vanishing a and increasing Hc2{0) with increasing disorder [|,|l§ 

In the case of Pt-substitutions in the investigated range < y < 0.75, the values of Tc, (3 
and Hc2{0) are reduced, too. As for R-substitutions, those superconducting properties exhibit 
minima at intermediate composition while the Sommerfeld constant and the curvature 
parameter a of Hc2{T) depend monotonically on y. For y < 0.2 a strong decrease of a 
with increasing y is observed, but for y > 0.2 an increase of a does not occure (see fig. 
This behaviour of Hc2 suggests that the quasi-dirty limit has been reached at about y « 0.2 
since i?c2(0) increases above y « 0.2, while the curvature of Hc2 measured by a remains 
strongly reduced. The results obtained for Y(Nii„j,Ptj,)2B2C show that the deviations from 
the linearity of "f{H) measured by (3 are not correlated with the field exponent a. While a 
almost vanishes, /3 does increase for y > 0.2. Thus, here the behaviour of the specific heat in 
the vortex state even in the quasi-dirty limit remains rather complex. In this context a similar 
behaviour for (Nb,Ti) alloys which are apparently strongly disordered superconductors is 
noteworthy: on one hand, a strong sublinear j{H) oc i/^'^-dependence is observed, and on the 
other hand, a standard parabolic Werthammer-Helfand-Hohenberg temperature dependence 
of Hc2{T) occures, i.e. a negative curvature near Tc is present p5| . 

To summarize, we have shown that the "f{H) curves of the pure specimens {x = 0; 1) 
of Ya;Lui_a;Ni2B2C exhibit the strongest sublinear behaviour reported for superconductors. 
Weak disorder effects caused by isoelectronic substitutions of Lu by Y yield a reduction of 
the 7 (i?)-non- linearity without reaching the standard linear behaviour. Similar moderate 
suppressions of characteristic features which are typical for the quasi-clean limit have been 
found for the upper critical field iJc2(0), the curvature exponent a, ■^n and Tc. Stronger 
disorder effects are caused by isoelectronic substitutions of Ni by Pt. From the behaviour 
of Hc2{T) it is deduced that a transition from clean to quasi-dirty limit occures caused by 
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0.0 0.5 1.0 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Y-concentration x Pt-concentration y 



Fig. 4 — Composition dependence of the transition temperature Tc derived from the onset of the jump of 
Cp(T) (upper panels), the fiommerfeld constant 7]v (second panels), the upper bound for the upper critical 
field H*2 according to eq. (W) and -ffc2(0) according to eq. jH) (third panels, see text for more details), the 
specific heat curvature exponent /3 of 7(H) according to eq. (H (fourth panels), and the curvature exponent a 
of the upper critical field Hc2 according to eq. (kh (lower panels) determined for Y2!Lui_2:Ni2B2C (left ) and 
for Y(Nii_j,Ptj,)2B2C (right). The lines are guides to the eye. 



isoelectronical substitutions at the T-site. The quasi-dirty hmit is deduced from the nearly 
vanishing curvature of Hc2{T) for y > 0.2. At the same time there the subhnearity of ^{H) 
remains and does even increase. Hence, a simple monotonical relationship between a and 
P, as one might expect by considering the results on Y2,Lui_2;Ni2B2C only, does not hold 
in the quasi-dirty limit. In the case of intermediate deviations from the linearity of ^{H) 
{(3 = 0.2 — 0.35) our results on specific heat at low magnetic fields are discussed in the context 
of a dirty d-wave model on the one hand and within the framework of the conventional s-wave 
picture in the quasi-clean limit on the other hand. At low fields the H\nH dependence of 'y{H) 
predicted for d-wave pairing in the dirty (unitary) limit is not very distinct from the H^~^ 
behaviour which favoures s-wave superconductivity in the quasi-clean limit. Thus, results on 
7(i7) only are insufficient to rule out unconventional pairing in borocarbide superconductors. 
But from the few available Cp data points at low-temperature and i? = a small gap of 
A(0)/A;bTc w 0.5 might be derived for a clean YNi2B2C sample. This reveals additional 
support for an anisotropic s-wave pairing state in non-magnetic borocarbides. 



* * * 
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